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^ . Abstract 
^vq ' Large discrepancies between quasi-free neutron-neutron (nn) cross section data from neutron-deuteron 

(nd) breakup and theoretical predictions based on standard nucleon-nucleon (NN) and three-nucleon (3N) 
forces are pointed out. The nn Sq interaction is shown to be dominant in that configuration and has to be 
increased to bring theory and data into agreement. Using the next-to- leading order (NLO) Sq interaction 
of chiral perturbation theory (%PT) we demonstrate that the nn QFS cross section only slightly depends 
on changes of the nn scattering length but is very sensitive to variations of the effective range parameter. 
In order to account for the reported discrepancies one must decrease the nn effective range parameter by 
about w 12% from its value implied by charge symmetry and charge independence of nuclear forces. 

1 PACS numbers: 21.45.-v, 21.45.Bc, 25.10.+S, 25.40. Cm 
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I. INTRODUCTION 



The knowledge of the NN interaction is fundamental for interpreting nuclear phenomena. The 
proton-proton (pp) experiments provide a solid data basis [l[ [J, which restricts theoretical as- 
sumptions about the strong part of the pp force. In case of the neutron-proton (np) system this is 
only true to a smaller extent. The partial wave analysis of the np data [2] relies on the assumption 
that the isospin t = 1 piece can be taken over from the pp system and only the t = part is 
free in the adjustment to the data. The lack of a free neutron target forbids neutron- neutron (nn) 
experiments, therefore the information on the nn interaction can be deduced only in an indirect 
way. To that aim the best tool seems to be the study of the three-nucleon (3N) system composed 
from two neutrons and the proton. It is simple enough to allow a rigorous theoretical treatment, 
e.g. in the framework of Faddeev equations 3J. The neutron-deuteron elastic scattering together 
with the neutron induced deuteron breakup, supplemented with the triton properties, offer a data 
basis which can be used to test properties of the nn force. Especially the nd breakup process with 
its rich set of configurations for free three outgoing nucleons seems to be a powerfull tool to test 
the nuclear Hamiltonian. By comparing theoretical predicitons to the nd breakup data in different 
configurations not only can the present day models of two-nucleon (2N) interactions be tested, but 
also effects of three-nucleon forces (3NF's) can be studied. 

The nn quasi-free scattering (QFS) refers to a situation where the outgoing proton is at rest in 
the laboratory system. In the nd breakup also np QFS is possible. Here one of the neutrons is at 
rest while the second neutron together with the proton form a quasi-freely scattered pair. 

The reported nn QFS cross sections taken at E l £ b = 26 MeV [4] and at E l £ b = 25 MeV 
overestimate the nd theory by ~ 18%. Surprisingly, when instead of the nn pair the np pair 
is quasi-freely scattered, the theory follows nicely the np QFS cross section data taken in the 
E l £ h = 26 MeV nd breakup measurement [4]. That good description of the np QFS cross section 
contrasts with the drastic discrepancy between the theory and the nn QFS cross section data taken 
in the same experiment 14 1. 

nn 

We do not expect surprises in case of the pp QFS data 16H8J, since the information of the rich 

set of pp data has been incorporated into the pp forces. In fact a recent analysis jij] including 

the Coulomb force to pp QFS data lead to a nice agreement, while in previous analysis f]-8] the 

Coulomb force was not yet included. Additional theoretical efforts to include all effects of the 

I I 

Coulomb force beyond the ones in [9|] are underway. 

In section |TT] we exemplify the stability of the QFS cross sections against changes of modern 
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nuclear forces. We also demonstrate that below ~ 30 MeV the 1 Sq and 3 Si- 3 -Di NN force compo- 
nents dominate the QFS cross sections. In section Hill we analyse the np as well as the nn QFS data 

n 

from |4] in terms of rigorous solutions of the 3N Faddeev equation and discuss necessary changes 
in the 1 Sq nn force component to remove the discrepancies in the nn QFS cross section. Thereby 
a detailed study is performed using the next-to-leading order (NLO) chiral NN force, composed 
of contact interactions and the one-pion exchange potential. It reveals that the effective range 
parameter is decisive to reconcile theory and data. The outcome is discussed in section [IV] and 
further experimental insights on the nn force are proposed. Finally we summarize in section IVl 



II. STABILITY AND SENSITIVITY STUDIES 



It is known that nd scattering theory provides QFS cross sections which are highly independent 



from the realistic NN potential used in the ca 



culations and that they practically do not change 



when any of the present day 3NF's is included 3|, 



id, 



111 ]. We exemplify it in Fig. [TJ for the nn and 



np QFS geometries of ref. [4|. There results of 3N Faddeev calculations [3J based on different high 
precision NN forces (CD Bonn Nijm I and Nijm II [lj]) alone or combined with the TM99 
3NF Q,Q are shown. 



The sensitivity study performed in [lOj revealed that at energies below ps 30 MeV the 1 S and 
3 S\- 3 Di NN force components provide the most dominant contribution to the QFS cross sections 
with much less contributions of higher partial waves. Specifically, in the np QFS geometries the 
3 S'i- 3 Di is the dominant force component while for nn QFS it is the 1 Sn force which contributes 
decisively. Again we exemplify it for nn and np QFS geometries of ref. [J] in Fig. [21 

Such a dominance for the QFS peak is understandable since the QFS cross sections are prac- 
tically insensitive to the action of the presently available 3NF. Then at low energy the largest 
contribution should be provided by the S-wave components of the NN potential. In case of free 
np and nn scattering these are the 1 5o(np)+ 3 S'i — 3 D\ and 1 S'o(nn) contributions, respectively. In 
the simple minded spirit that under QFS condition one of the three nucleons (at rest in the lab 
system) is just a spectator such a dominance of a two-nucleon encounter is to be expected. In 
reality, however, the projectile nucleon also interacts with that " spectator" particle and the three 
nucleons at low energies undergo higher order rescatterings [3, 12]. Thus the scattering to the 
final nn (np) QFS configuration also receives contributions from the np 3 Si — 3 D\ (nn 1 Sq) inter- 
action. Despite of all that the numerical results clearly reveal that for the np QFS configuration 
the 3 Si — 3 D\ force is the most dominant contribution and for the nn QFS it is the 1 S'o force (for 
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free nn scattering there is no 3 Si — 3 D\ interaction possible). This implies that the nn QFS is a 
powerful tool to study the 1 So nn force component. 

That extreme sensitivity of the nn QFS cross section to the 1 So nn force component is demon- 
strated in Fig. [3] for the QFS geometries of ref. (4]. To that aim we multiplied the 1 So nn matrix 
element of the CD Bonn potential by a factor A. The result is, that the nn QFS cross section 
undergoes significant variations while the np QFS cross section is practically unchanged. The dis- 
played A-parameters include also the value A = 1.08 which is necessary to get agreement with the 
nn QFS data of ref. El- 

While both, 1 So and 3 Si — 3 D\, np forces are well determined by np scattering data (with the 
restrictions mentioned above) and by the deuteron properties, the 1 So nn force is determined up 
to now only indirectly due to lack of free nn data. The disagreement between data and theory in 
the nn QFS peak points to the possibility of a flaw in the nn 1 So force. It was shown in 10] that 
in order to remove the ~ 18% discrepancy found in [4| for the nn QFS cross section required an 
increased strength of the 1 So nn interaction which when given in terms of a factor A amounts to 
A « 1.08. In Fig.H we show the effect of the A-modification for the nn scattering length CL nn and 
for the effective range parameter r e jf, and in Fig. [5] for the binding energy of two neutrons in the 
1 Sq state. It is seen that taking A = 1.08 leads to a nearly bound state of two neutrons. 



III. IMPLICATIONS ON THE ^0 NN EFFECTIVE RANGE PARAMETER 

Since the multiplication of the 1 Sq potential matrix element by a factor A induces changes in 
the effective range as well as in the scattering length the question arises, which from both effects is 
more important for the nn QFS cross section variations ? To answer that question we performed 
3N Faddeev calculations based on the next-to-leading (NLO) order xPT potential [ltIQ] including 
all np and nn forces up to the total angular momentum j ma x = 3 in the two-nucleon subsystem. 
The 1 Sq component of that interaction is composed of the one-pion exchange potential and contact 
interactions parametrized by two parameters Ci g and Ci So 

v( 1 s ) = a So + a So (p 2 +p 2 ) ■ (i) 

Standard values are Ci So = -0.1557374 * 10000 GeV" 2 and Ci So = 1.5075220 * 10000 GeV" 4 for 
cut-off combinations {A, A} = {450 MeV , 500 MeV} [IS]. 

Multiplying Cig by a factor C^So) and Cig by a factor Ci( 1 5o) one can induce changes of 
the nn 1 S'o interaction. Requiring either the value of the scattering length a nn or the value of the 
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effective range parameter r e fj to be constant correlates the Ci( 1 S'o) and C^So) factors. 

Changing Ci( 1 S'o) and C^So) in such a way that the scattering length is kept constant and 
equal a nn = —17.6 fm leads to changes of the effective range r e ff shown in Fig. [6] The resulting 
changes of the nn and np QFS cross sections for geometries of ref. [J] are shown in Fig. [7J for five 
sets of Ci( 1 S'o) and C^So) factors with different nn 1 Sq effective range parameters ranging from 
r e ff = 2.03 fm to r e fj = 3.07 fm; one of them corresponding to the value required by the data. 

Similarily, changing Ci( 1 5'o) and C^So) while keeping the effective range constant to r e fj = 
2.75 fm, leads to changes of the nn 1 Sq scattering length a nn shown in Fig. [8] The resulting changes 
of the nn and np QFS cross sections are presented in Fig. [9] for four values of the nn 1 Sq scattering 
length ranging from a nn = —10.9 fm to a nn = —75.9 fm. It is clearly seen that the nn QFS cross 
sections depend only slightly on a change of the scattering length. The variations of the QFS cross 
section maximum stays below ~ ±4%. On the other side much stronger variations of the nn QFS 
cross sections result from changes of the effective range (see Fig. [7]). 

Thus we can conclude that the A-enhancement mechanism for the 1 S'o nn force studied in [l 
acts mainly through the change of the effective range parameter. Thus in order to remove the 
discrepancies found in Q and for the nn QFS cross section a change of the nn So effective 
range parameter is required. Its value taken under the assumption of charge symmetry and charge 
independence of nuclear forces is r e ff = 2.75 fm and it has to be changed to r e jf ~ 2.41 fm. That 
implies a large charge symmetry and charge independence breaking effect of about ~ 12% for that 
parameter. 

We would like to add that the discussed changes of r e ff did not affect the elastic nd cross 
section nor vector or tensor analysing powers to a measurable extent. Only more complicated spin 
observables in elastic nd scattering are affected but the present day experimental errors are much 
larger than those changes. 

IV. DISCUSSION AND FURTHER EXPERIMENTAL INFORMATION 

Is such a large isospin breaking effect at all possible in view of the present understanding of 
nuclear forces? 

First of all it seems unprobable that only the effective range would reveal large isospin breaking 
and the scattering length will be left unaffected. In xPT the leading isospin breaking contribution 



is provided by isospin breaking contact interaction without derivatives 



191 ] . It turns out that the 



effective range parameter is quite insensitive to that isospin breaking contact force and typical 
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isospin breaking effects for r e ff are small and under ~ 1% [191 ] . 

The reported discrepancies for nn QFS require however a much larger effect for r e ff of the order 
~ 12%. Only when the contact terms in next orders would be unnaturally large one could expect 
larger isospin breaking effects for r e ff. Assuming naturalness it seems rather unprobable. 

Since it seems unlikely that isospin breaking effects will show up, if at all, in the effective 
range parameter alone without affecting simultaneously the nn scattering length, the question of 
a possible existence of a bound state of two neutrons reappears. 

Present day NN interactions allow only one bound state of two nucleons, namely the deuteron, 
where the neutron and the proton are interacting in a state with angular momenta I = or 2, total 
spin s = 1, and total angular momentum j = 1. When the neutron and proton are interacting 
with the 1 S'o force no bound state exists and only a virtual resonant state occurs as documented 
by the negative scattering length a np = —21.73 fm. Also the data for the proton-proton system 
exclude a 1 Sq pp bound state; however in this case the nuclear force is overpowered by the strong 
pp Coulomb repulsion. Assuming charge-independence and charge-symmetry of strong interactions 
also the two neutrons should not bind in the 1 So state. 

It also seems that modern nuclear forces do not allow for the 3n and 4n systems to be bound 
[20! . However, in view of the strong discrepancies between theory and data found in the nd breakup 
measurements for the nn QFS geometry, which cannot be explained by present day nuclear forces, 
it appears reasonable to check experimentally the possibility of two neutrons being bound. 

There are reactions which provide conditions advantegous for a hypothetical di-neutron bound 
state. Such conditions can be found e.g. when two neutrons are moving with equal momenta 
and with relative energy close to zero. That occurs in the so called final-state-interaction (FSI) 
geometry of the nd breakup. Incomplete nd breakup measurements have been performed in the 
past to study properties of the 1 5q nn force 21[. Even a dedicated experiment was performed in 



22| in which the spectrum of the proton going 



order to look for a hypothetical 1 Sq nn bound state 
in forward direction has been measured with the aim of a precise determination of its high energy 
region. The negative result of {22] showed that the nd reaction is not suitable for such a study. 

It seems that much more appropriate would be reactions in which from the begining two neutrons 
occupy a configuration advantegous for their binding. 

It is known {23]) [2^ that 3 He is predominantely a spatially symmetric S state with its two 
protons mainly in opposite spin states. This component amounts for ~ 90% of the 3 He wave 
function. Similarily, the two neutrons in 3 H are restricted to be in a spin-singlet state. That makes 
the triton target a very suitable tool to look for a nn bound state in 7 induced breakup of 3 H. The 
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idea is to measure the spectra of the outgoing protons in such a reaction. The two-neutron bound 
state, if existant, should reveal itself as a peak above the highest available proton energy from the 
3-body decay of 3 H. We show in Figs. [10] to [11] the outgoing proton spectra from the , y( 3 H,p)nn 
reaction for a number of 7 energies and angles of the outgoing protons. The big advantage of that 
reaction is that 7 interacts predominantely with the proton. 

Also other reactions, such as e.g. 3 H(n,d)nn and 3 H(d, 3 He)nn, provide conditions advantegous 
for two neutrons to bind. They are complimentary and independent from the 3 H(7,p)nn reaction 
and the data from all three processes should provide an answer to the question whether two neutrons 
can form a bound state. The reaction 3 H(d, 3 He)nn cannot presently be treated in a theoretically 
rigorous manner, however with the rapid increase in computer power such a treatment based on 
Fadeev-Yakubovsky equations can be expected in the near future. 



V. SUMMARY 

The strong discrepancy in the nn QFS nd break up configuration found in [J, is reconsidered. 
It is documented again that at low energies (below ~ 30 MeV) the nn (np) QFS cross section 
depends dominantly only on the l So ( 3 Si — 3 D\) NN force component and higher partial wave 
contributions are quite small. Furthermore the theoretical results are quite stable under exchange 
of the standard nuclear forces. Also the present day available 3N forces have negligible effect on 
the QFS configurations. Since no direct measurement of the nn force is available there is the 
possibility that the properties of the nn force are still unsettled. Thus simply multiplying the nn 
1 5o force matrix element by a factor A = 1.08 one can perfectly well reconcile theory and data. In 
addition we performed a more detailed study using the NLO chiral potential, which is composed 
of the one-pion exchange and contact interactions depending on two parameters. That dependence 
allowed us to study separately variations in the scattering length a nn leaving the effective range 
parameter r e ff constant and vice versa. Thereby it turned out that the nn QFS peak height is 
very sensitive to r e fj and hardly sensitive to a nn . The outcome for an agreement with the data 
is the requirement that r e ff decreases from the value r e ff = 2.75 fm to a significantly smaller 
one, r e jf = 2.41 fm. That strongly breaks charge symmetry and charge independence and is not 
supported by present day chiral potential theory. 

So, what might be a solution to remove the discrepancy? 

If the data are taken for granted there remains the possibility that a di-neutron exists. We pro- 
pose additional experimental investigations, like the 3 H('y,p)nn process and evaluated the proton 
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spectra at various emission angles emphasizing its high energy region. 

The direct inclusion of A-degrees of freedom into allows for a rich set of additional NN 
and 3N force diagrams which are presently under investigation [2^]. This might reconcile theory 
and data also for the space-star discrepancy |3J in the nd breakup process. 

Right now the situation is unsettled. 
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FIG. 1: (color online) The cross section d 5 a / dflid^dS for the E\f h — 26 MeV nd breakup reaction d(n, nn)p 
(upper panel) and d(n, np)n (lower panel) as a function of the S-curve length for two complete configurations 
of Ref. 4]. QFS nn refers to the angles of the two neutrons: Q\ = 02 = 42° and QFS np refers to the angle 
B\ = 39° of the detected neutron and 82 = 42° for the proton. In both cases (£12 = 180°. The (practically 
overlapping) lines correspond to different underlying dynamics: CD Bonn [13] - dashed (blue), Nijm I - 
dotted (black), Nijm II [14] - dashed-dotted (green), CD Bonn+TM99 - solid (red), Nijm I +TM99 Q,Q 
- dashed-double-dotted (orange), Nijm II + TM99 - double-dashed-dotted (maroon). All partial waves with 
2N total angular momenta up to j ma x = 5 have been included. 
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FIG. 2: (color online) The cross section d 5 a / ' dQid^dS for the E l £ h — 26 MeV nd breakup reaction d(n, nn)p 
(upper panel) and d(n, np)n (lower panel) as a function of the S-curve length for two complete configurations 
of Ref. 0| specified in Fig. [TJ The different lines show contributions from different NN force components. 
The solid (red) line is the full result based on the CD Bonn potential [13J and all partial waves with 2N total 
angular momenta up to j ma x = 5 included. The dotted (black), dashed-dotted (green), and dashed (blue) 
lines result when only contributions from 1 S'o, 3 5i — 3 D\, and 1 S'o + 3 Si — 3 D\ are kept calculating the 
cross sections. The dashed-double-dotted (brown) line presents the contribution of all partial waves with 
the exception of 1 S$ and 3 <Si — 3 D\. 
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FIG. 3: (color online) The cross section d 5 a / dflid^dS for the E l ^ h — 26 MeV nd breakup reaction d(n, nn)p 
(upper panel) and d(n, np)n (lower panel) as a function of the S-curve length for two complete configurations 
of Ref. specified in Fig. [TJ The lines show sensitivity of the QFS cross sections to the changes of the 
nn Sq force component. Those changes were induced by multiplying the 1 Sq nn matrix element of the 
CD Bonn potential by a factor A. The solid (red) line is the full result based on the original CD Bonn 
potential [13| (a nn = —18.8 fm, r e ff = 2.79 fm) and all partial waves with 2N total angular momenta up 
to j m ax — 5 included. The dashed (blue), dotted (black), and dashed-dotted (green) lines correspond to 
A = 0.9 (a„„ = -8.3 fm, r eff = 3.12 fm), 0.95 (a„„ = -11.7 fm, r eff = 2.96 fm), and 1.05 (a„„ = -42.0 fm, 
r e ff — 2.66 fm), respectively. The double-dashed-dotted (violet) line shows cross sections obtained with 
A = 1.08 (a n n = —134.7 fm, r e // — 2.61 fm), which factor is required to get agreement with nn QFS data 
of ref. Ul. 
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FIG. 4: (color online) The changes of the nn scattering length a nn and the effective range parameter r e ff 
with factor A by which the So nn matrix element of the CD Bonn potential is multiplied: V r „„( 1 S'o) = 

A * Vc'D Bonn^So). 




FIG. 5: (color online) The range of A values by which the 1 Sq nn matrix element of the CD Bonn potential is 
multiplied (V nn ( So) = \*Vcd Bonn{ 1 So)), for which the two neutrons form a bound state with the binding 
energy E b . 
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FIG. 6: (color online) Changes of the effective range parameter r e // in the 1 Sq partial wave caused by 
a correlated change of the factors Ci( 1 5'o) and C^So) as shown in bottom part of this figure. This 
correlation between the factors Ci( 1 S'o) and C^So) corresponds to a constant value of the scattering 
length a nn = —17.6 fm. 
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FIG. 7: (color online) Changes of QFS cross sections for configurations specified in Fig.[T]caused by correlated 
change of factors Ci( 1 5o) and C^So) shown in Fig. [6] All lines show results of Faddeev calculations based 
on LO xPT potential and all partial waves with 2N total angular momenta up to j max — 3 included. 
They differ in the nn So force which was obtained keeping constant scattering length a nn = —17.6 fm 
and changing constants Ci( 1 S'o) and C2( 1 S'o) to get different effective ranges which are: solid (red line) 
- CxC-Sq) = 1.0, C 2 ( 1 S ) = 1.0, r eff = 2.75 fm, dashed (blue line) - C^Sq) = 1.5, C 2 ( 1 S ) = 1.0615, 
r eff = 3.07 fm, dotted (black line) - C^Sq) = 0.8, C 2 { 1 S ) = 0.9275, r eff = 2.54 fm, dashed-dotted (green 
line) - Ci^So) = 0.5, C 2 CS ) = 0.7675, r eff = 2.03 fm. The double-dashed-dotted (violet) line shows 
cross sections obtained with Ci( 1 S'o) = 0.7064, C , 2( 1 S'o) = 0.8842, r e ff = 2.41 fm, which are required to get 
agreement with nn QFS data of ref.Q. 



15 



d 
d 



•100 



c3 



-200 

1,2 



1,0- 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 


0,9 1,0 1,1 

c 2 ('s o) 



0,9 



1,2 

c 1 (X) 



1,5 



FIG. 8: (color online) Changes of the nn scattering length a nn in the 1 5o partial wave caused by a correlated 
change of the factors C\(}Sq) and C2CS0) as shown in the bottom part of this figure. This correlation 
between the factors Ci( 1 S'o) and C , 2( 1 <S'o) corresponds to a constant value of the effective range parameter 
r e ff = 2.75 fm. 
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FIG. 9: (color online) Changes of QFS cross sections for configurations specified in Fig. Q] caused by a 
correlated change of the factors Ci( 1 S'o) and C , 2( 1 5o) shown in Fig. [8j All lines show results of Faddeev 
calculations based on the NLO xPT potential and all partial waves with 2N total angular momenta up to 
jmax = 3 included. They differ in the nn 1 Sq force which was obtained keeping the effective range parameter 
r e ff = 2.75 fm constant and changing the constants Ci( 1 S'o) and C^So) to get different scattering lengths 
which are: solid (red line) - C^Sq) = 1.0, C 2 ( 1 S ) = 1.0, a nn = -17.6 fm, dashed (blue line) - Ci^So) = 
0.8, C2C-S0) = 0.8953, a nn = -10.9 fm, dotted (black line) - Ci^Sq) = 1.3, C 2 ( 1 S ) = 1.1139, a nn = 
-45.3 fm, dashed-dotted (green line) - C^Sq) = 1.4, C 2 ( 1 S Q ) = 1.1410, a nn = -76.0 fm. 
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FIG. 10: (color online) The spectra of the outgoing proton from the reaction 3 H( r y, p)nn with E 1 = 10 MeV 



at different lab. angles of the proton. They have been calculated using the AV18 



the current composed of single nucleon and meson exchange currents 
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FIG. 11: (color online) The spectra of the outgoing proton from the reaction 3 i? (7, p)nn with E 1 = 10 MeV 



at different lab. angles of the proton. They have been calculated using the AV18 



the current composed of single nucleon and meson exchange currents 
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